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Abstract 

We study charge excitation spectra in the two-dimensional t-J model on a square lattice to explore a charge-order 
tendency recently found in electron-doped cuprates around the carrier density 0.15. The static susceptibility of d- 
wave charge density, which corresponds to the nematic susceptibility at the momentum transfer q = (0, 0), shows two 
characteristic peaks at momenta of the form qi = and q 2 = (9,0). These two peaks originate from the so- 

called 2kF scattering processes enhanced by the d-wave character of the bond-charge density. The peak at qi is much 
broader, but develops to be very sharp in the vicinity of its instability, whereas the peak at q 2 becomes sharper with 
decreasing temperature, but does not diverge. The equal-time correlation function, which is measured by resonant 
x-ray scattering, exhibits a momentum dependence similar to the static susceptibility. We also present energy-resolved 
charge excitation spectra. The spectra show a V-shaped structure around q = (0, 0) and bend back toward close to 
zero energy due to the charge-order tendency at qi and q 2 . The resulting spectra form gap-like features with a 
maximal gap at q ~ qi/2 and q2/2. We discuss implications for the recent experiments in electron-doped cuprates. 


PACS: 74.72.Ek Electron-doped cuprates', 

75.25.Dk Orbital, charge, and other orders, including coupling of these orders', 
78.70.Ck X-ray scattering 


1 Introduction 

Charge order (CO) in high-temperature cuprate supercon¬ 
ductors attracts renewed interest. CO is known in La- 
based materials as a spin-charge stripe order [I] , in which 
CO is accompanied by a spin order. However, a different 
type of CO has been observed recently in various hole- 
doped cuprates such as Y- EnmiBiaiziiH], Bi- [iiiniiii], 
and Hg-based [12] materials. In these materials, the CO is 
not accompanied by a spin order. Furthermore a modula¬ 
tion vector of the CO decreases with doping, the opposite 
tendency observed in the La-based materials. The origin 
of the newly found CO as well as its relation to supercon¬ 
ductivity is under active debate naiTiiiniiiiKn]. 

Research interest also goes to electron-doped cuprates. 
Resonant inelastic x-ray scattering (RIXS) reveals that 
charge excitation spectra develop to form a V-shaped dis¬ 
persion m around the momentum q = (0, 0) and extends 
up to around 1.5 eV at q = (O.Gtt, 0) and (O.Gtt, O.Btt) [12] 
in Nd 2 -a;Cea;Cu 04 with x = 0.15. Quite recently resonant 
x-ray scattering (RXS), which integrates a RIXS spectrum 


up to infinity with respect to energy, has revealed a charge 
excitation peak at q ~ (0.487r, 0) near a; = 0.15 [2D]. The 
observed wavevector is rather close to that found in hole- 
doped cuprates [21[l[l[51[i[71IHl|Dl[ini[IIl[Il], implying 
a possible universal phenomenon for the CO in cuprate 
superconductors. The correlation length of the CO is, 
however, estimated to be 4-7 lattice constants, i.e., it is 
not a long-range order [2D]. 

The wavevector of q ss (0.487r, 0) obtained by RXS [2D] 
is covered by the RIXS by Ishii et al. [ID] , but the observed 
RIXS spectra do not seem to suggest clearly some char¬ 
acteristic feature associated with a CO, such as softening 
of the spectrum toward a long-range order at the corre¬ 
sponding wavevector. This peculiar situation motivates 
us to study charge excitations in electron-doped cuprates 
more closely from a theoretical point of view. 

Charge excitations in electron-doped cuprates are not 
much known theoretically. Ishii et al. studied the usual 
density-density correlation functions and discussed the 
RIXS spectra m- Bejas et al. EH, on the other hand, 
studied all possible COs in the t-J model with parameters 
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appropriate for electron-doped cuprates. They found that 
instead of a usual charge-order instability, various types 
of bond order tend to occur much more strongly. In par¬ 
ticular, a d-wave bond-order tendency is dominant in a 
moderate doping region. While its instability is expected 
at a wavevector close to (tt, tt) , they found a meta-stable 
solution of the d-wave bond order at q Ri (0.497r, 0). This 
wavevector is very close to the experimental observation 
by RXS [10]. 

Encouraged by this agreement with the experiment, 
we study charge excitations associated with a d-wave 
bond order in the two-dimensional t-J model on a square 
lattice by taking parameters appropriate for electron- 
doped cuprates. We compute three quantities: static d- 
wave bond-order susceptibility Xd(q, 0), its spectral weight 
Imyd(q, w), and the corresponding equal-time correlation 
function S{q). The second and third quantities can be 
measured directly by RIXS and RXS, respectively. Our 
obtained results capture essential features observed in ex¬ 
periments such as a V-shape dispersion of Imxd(q, w) 
near q = (0,0) [111 HH and a short-range CO with 
q r: (0.487r,0) [20]. In addition, we obtain several new 
insights: First, a CO is expected also at qi = {q',q') with 
q' R 0.847r. In fact this CO has a stronger intensity than 
the CO at q 2 = (g, 0) with q r 0.497r. However, the peak 
at qi is much broader in momentum space than that at 
q 2 and becomes sharp only in the vicinity of its instabil¬ 
ity. Second, the dispersive peak of Imyd(q, w) bends back 
toward close to zero energy at qi and q 2 where Xd(q, 0) 
and S'(q) exhibit a peak. The resulting charge excitation 
spectra show gap-like features between qi and q = (0, 0), 
and between q = (0,0) and q 2 , with a maximal gap at 
q R iqi and dq2. 


2 Model and formalism 


Various approximations to the t-J [H] ES] [21] and the 
strong coupling Hubbard [24] model show that the mod¬ 
els have a strong tendency toward phase separation, es¬ 
pecially for parameters appropriate for electron-doped 
cuprates. The phase separation, however, can be an ar¬ 
tifact caused by neglecting the long-range Coulomb in¬ 
teraction. In fact, the Coulomb interaction term appears 
naturally when the t-J model is derived from the three- 
band Hubbard model [IS]. Hence we include the nearest- 
neighbor Coulomb interaction in the t-J model as a min¬ 
imal model to study electron-doped cuprates. Our model 
then reads 


H — 


cLc, 


■ + jY^ 

do) 


Si * Sj ^TliTlj 


(iJ) 


( 1 ) 

where tij = t {t') is the hopping between the first (second) 
nearest-neighbor sites on a square lattice, J and V are the 
exchange and Coulomb interactions between the nearest- 
neighbor sites, respectively. (i,j) indicates a nearest- 


neighbor pair of sites. c|^ {cia) is the creation (annihi¬ 
lation) operator of electrons with spin a {a =4';t) in the 
Fock space without any double occupancy, ni = 
is the electron density operator and Si is the spin opera¬ 
tor. The role of the V-term turns out to merely suppress 
phase separation in a doping and temperature (T) range 
relevant to cuprates. In fact, our obtained results are not 
affected essentially by the V-term. 

In leading order of a 1 /A^-expansion [26] . the kinetic 
term of the electrons is characterized by an effective elec¬ 
tronic dispersion 


Ek = —2 


t^ + JA 


( 2 ) 

where 6 is the doping rate away from half-filling and /i is 
the chemical potential. The bare hopping integrals t and t' 
are renormalized by a factor of S. The term JA in eq. @, 
which is not present at bare level, comes from the exchange 
term (the second term in eq. ([IJ). The magnitude of A 
describes a bond amplitude between the nearest neighbor 
sites. Values of A and ^ are determined self-consistently 
at a given 6 by solving the following equations: 


A = 


4A, 


(cos kx + cos ky)np(eii,) 


(3) 


and 

l-<^= (4) 

k 


Here riF is the Fermi function and Ns the total number of 
lattice sites. 

In the above scheme Bejas et al. studied all possible 
charge instabilities for both hole-doped[T7] and electron- 
doped[3I] cuprates. They found that a relevant instability 
to discuss the electron-doped cuprates around 5 = 0.15 is 
a d-wave bond order where bond amplitude is modulated 
along both x and y direction, and its relative phase is in 
antiphase 0. This ordering pattern is shown in fig. [T] by 
choosing wavevectors close to those relevant in the present 
study. 

In the present study, we explore closely a charge-order 
tendency associated with the d-wave bond order. In par¬ 
ticular, we aim to give some insight into the charge exci¬ 
tations recently observed by RIXS and RXS in electron- 
doped cuprates from a theoretical point of view. Following 
m and [21] we focus on the effective dynamical d-wave 
charge susceptibility. 


Xd(q,w) 


(8JA^)-^ 

1 — 2jn(q, w) 


(5) 


which becomes exact in leading order of 1/A. The bare 
polarizability n(q, w) reads 


n(q,a;) 


1 2|'u\ ^^(^k+q/ 2 ) ^F(ek-q/2) 

^ ^ ek+q/2 - ek-q/2 - UJ - 11] ’ 


^This order is refereed to as BOPa:y in m 
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Figure 1: (Color online) d-wave bond order for qi w q = 
(O.Stt, O.Stt) (a) and q 2 ~ q = (0.57r,0) (b). The black 
lines denote a stronger (solid line) and weaker (dotted line) 
bond relative to the average bond amplitude (gray line), 
which corresponds to A in eq. ()3|). The width of the lines 
indicates the modulation amplitude, (c) 2kp scattering 
processes, which determine the wavevectors qi and q 2 . 

where ry(> 0) is an infinitesimally small value and we take 
rj = 10“^ because of a practical reason of numerical com¬ 
putations. The form factor 7 (k) = {cosk^ — cosfcy)/2 
comes from the intra-unit-cell symmetry. This form fac¬ 
tor has d-wave symmetry, and corresponds to the 

well known nematic susceptibility for q = 0 [271 EH [29]. 
The property of Xd(q, w) near q = 0 was already studied in 
[30] by focusing on a collective mode of the d-wave bond 
order in both paramagnetic and superconducting states. 
Here we study eq. in a different situation in which 
COs tend to occur at q = qi and q 2 . 

In what follows we present results for J/t = 0.3 and 
t'/t = 0.30, which are suitable for electron-doped cuprates. 
We fix the carrier density S = 0.15 so that our results will 
be compared directly with recent experiments [niEniEo]. 
Our conclusions do not depend sensitively on a precise 
choice of parameters and we choose V/t = 1. Below we 
present all quantity of the dimension of energy in units of 
t. 

3 Results 

We present the static susceptibility Xd(q, 0), the equal¬ 
time correlations function S'(q), and the spectral weight 
Imxd(q, w) in (q, w) space. RXS and RIXS measure S'(q) 
and Imxd(q, w), respectively. 

We first computed Xd(q, 0) in the entire Brillouin zone 
and found two well-defined peaks at qi = (0.847r, 0.847r) 



Figure 2: (Color online) The static susceptibility Xd(q, 0) 
along the directions of (0,0)-(7r, tt) (a) and (0, 0)-(7r, 0) (b) 
for various temperatures. In (a), the curve at T = 10“® is 
scaled by a factor of 0.5. 

and q 2 = (0.497r, 0) near zero temperature. To clarify their 
temperature dependence, we plot Xd(q, 0) along the (0,0)- 
(7r,7r) direction in fig. [2] (a). At T = 10“^ a very sharp 
peak forms at q = qi. This peak is due to the proximity to 
a quantum critical point of the d-wave bond-order insta¬ 
bility, which is present at Sc Ri 0.13[1I]. However, once the 
temperature is increased, the peak is immediately broad¬ 
ened and becomes less clear already at T — 0.01. In fig. 
|2](b) we plot Xd(q;0) along the (0,0)-(7r, 0) direction. In 
contrast to fig. [2] (a), the temperature dependence of the 
peak features more usual behavior in the sense that the 
peak is broad at high T and smoothly grows to be a pro¬ 
nounced peak at low temperature. In spite of this clear 
peak structure, X!i(qi0) does not diverges at q 2 . That is, 
there is no indication that the d-wave bond order becomes 
long range along the (0, 0)-(7r,0) direction. 

The static susceptibility Xd(q, 0) is a useful quantity 
to study the stability of a system, i.e., an ordering phe¬ 
nomenon. This quantity is, however, not measured di¬ 
rectly by RXS. Rather, RXS measures the equal-time cor¬ 
relation function <5'(q), which is defined by 

1 

5'(q) = - / dwlmxd(q,w) [?^b(w) + 1], (7) 

d-oo 

where np is the Bose factor. In fig. [3| (a) we show an 
intensity map of S'(q) along the (0,0)-(7r, tt) direction in 
a temperature range 0 < T < 0.1. Although the spectral 
weight tends to accumulate around q = qi with decreasing 
T, the temperature dependence is weak and the spectral 
weight still spreads down to zero temperature in spite of 
the proximity of the corresponding charge instability. To 
show the temperature dependence of S'(q) more clearly, 
we plot a spectrum AS'(q;T) = S'(q;T) — S'(q;T = 0.1) 
in fig. O (b). Its temperature dependence is very similar 
to that of the static susceptibility shown in fig. [2| (a) 
except at T = 10“^. In fig. |3| (c) and (d), we present 
the corresponding results along the (0,0)-(7r,0) direction. 
Although the CO tendency is stronger at q = qi than at 
q = q 2 (see fig. E]), the peak structure at q = q 2 is much 
clearer, being sharp and pronounced with decreasing T, 
as demonstrated in fig. |3|(d). 
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Figure 3: (Color online) (a) Intensity map of S'(q) along 
the (0, 0)-(7r, tt) direction for 0 < T < 0.1. (b) Evolution 
of the spectral weight AS'(q; T) = 5'(q; T) — S'(q; T = 0.1) 
for various temperatures. Consequently AS'(q;T) = 0 at 
T = 0.1. (c) and (d): corresponding results to (a) and (b), 
respectively, along the (0,0)-(7r, 0) direction. 


The energy-resolved spectral weight Imxd(q, w) is 
shown in fig. |4] at low temperature in the plane of q and 
oj. A V-shape dispersion develops from q = (0,0). This 
dispersion originates from individual particle-hole excita¬ 
tions and extends up to high energy. However, the spec¬ 
trum bends back and softens toward close to zero energy 
at q = qi and q 2 , where both static susceptibility (fig. 
El and equal-time correlation function (fig. [S]) exhibit 
a peak. These dispersion near qi and q 2 is interpreted 
as coming from collective charge excitations. This collec¬ 
tive feature is particularly clear near q = qi due to the 
proximity to the corresponding charge instability. A gap¬ 
like feature of charge excitations is visible between qi and 
q = (0, 0), and also between q = (0, 0) and q 2 , and forms 
a maximal gap of about 0.1 at q ~ ^qi and iq 2 . This 
gap-like feature is more pronounced along the (0, 0)-(7r, tt) 
direction because the d-wave character of the bond-charge 
density suppresses its low-energy scattering processes sub¬ 
stantially. 


4 Discussions 

Now we discuss implications for the experiments by 
RXSEO] and RIXS [^[H]. 

In fig. E] (c) and (d), the charge peak at q = q 2 be¬ 
comes sharper with decreasing temperature in the equal¬ 
time correlation function, but the real part of the suscepti¬ 
bility remains finite at the corresponding wavevector (fig. 
El (b)). We therefore conclude that the CO at q = q 2 
remains a short range, which is consistent with the RXS 



Figure 4: (Color online) Energy-resolved spectral weight 
Imy(i(q, w) in the plane of q and lo along the symmetry 
axes at low temperature. The spectral weight is scaled by 
a factor of 2 along the (0, 0)-(7r, 0) direction to get a better 
contrast to that along the (tt, 7r)-(0, 0) direction. 


measurements [20] . In particular, a short-range feature of 
the observed CO can be interpreted as an intrinsic prop¬ 
erty, i.e., it does not come from some disorders frequently 
present in actual materials. 

On the other hand, the present theory predicts that the 
CO tendency at q = qi is much stronger than at q = q2. 
This is because a long-range order with modulation vector 
q w qi occurs at T = 0 at the critical doping 5c ~ 0.13 
[21]. This peak at qi is, however, peculiar in the sense that 
it is much broader than at q = q2 and becomes sharper 
only in the vicinity of the charge instability. Since qi 
is rather close to (tt,tt), it is not straightforward to test 
it in experiment. In fact, it is difficult to perform RXS 
and RIXS up to near (7r,7r). Hence a usual x-ray diffrac¬ 
tion measurement can be more fruitful by exploring lattice 
modulations generated by the underlaying CO. Given that 
the CO instability at q ss qi is expected below 5 < ^cEIl, 
it might seem easier to measure a sample with lower car¬ 
rier density than 0.15. However, antiferromagnetism tends 
to be stabilized at lower carrier density and could mask 
the CO instability. 

The peak positions at qi and q 2 found in Xd(q, 0) (fig. 
Ej) and 5'(q) (fig. [3]) are determined mainly by two factors: 
the so-called 2kp scattering processes |3T] and the d-wave 
character of the bond order. The corresponding scattering 
processes are depicted in fig. [T] (c). In particular, the 
peak at q 2 becomes pronounced substantially by the d- 
wave form factor. Hence we expect that the observed CO 
at q = q 2 [20| has a d-wave character, which may be tested 
by RXS in the future [21]. 

As shown in hg. E] charge excitations feature a V- 
shaped spectrum around q = (0,0), which agrees quali¬ 
tatively with the experimental observations [Ulllg]. The 
spectra near q = (0,0) come mainly from individual 
particle-hole excitations, in favor of the experimental in¬ 
terpretation by Ishii et al. [19] . A quantitative compar¬ 
ison with the experiment requires additional care. The 
V-shaped dispersion reported in [19] extends to 1.5 eV 
at q = (0.67r,0) and (O.Ott, O.Gtt). Using t/2 = 500 meV 
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0, which is the estimated value for cuprates [33], our ob¬ 
tained dispersion (fig. |3|) extends up to w 500 meV at 
the same momenta. This energy scale is about a factor of 
three lower than the experimental observation. Our small 
energy scale originates mainly from a relatively small band 
width due to the renormalization of the bare t to an ef¬ 
fective hopping teff = as seen in eq. ([3|). On the other 
hand, a large band width develops immediately after dop¬ 
ing the Mott insulator phase, which cannot be captured 
quantitatively in terms of teff- In this sense, a quanti¬ 
tative comparison of energy scale of charge excitations is 
connected with a fundamental issue of doped Mott insu¬ 
lators and remains to be studied. Furthermore in present 
study we focus on charge excitations of d-wave bond order 
because it gives the most relevant contributions to charge 
excitations at low energy. For a full comparison with RIXS 
data, however, charge excitations from other types of bond 
orders m as well as the usual charge density should be 
considered since RIXS may contain spectra of those exci¬ 
tations especially in a high energy region. At low energy, 
contributions from the d-wave bond order should become 
dominant and thus it is interesting to test a softening of 
charge excitation spectrum at qi and q 2 (see fig. Sj) by 
RIXS and to clarify the actual energy scale there. 

The square lattice in the present model describes the Cu 
sites in the Cu 02 plane of cuprate superconductors and 
the center of the nearest-neighbor sites corresponds to the 
oxygen site. Therefore our bond order may be interpreted 
as a charge modulation at the oxygen sites [33] . 

It is natural to consider whether the present theory 
can be applied also to hole-doped cuprates by taking 
appropriate parameters. Comprehensive calculations in 
the hole-doped case m, however, did not capture a 
CO tendency compatible to the experimental observation 
[iiiiiiaisiiTiiiiiiiiniiiiKii]. Such calculations were 
performed in the paramagnetic state whereas in reality 
the CO is observed as an instability in the pseudogap 
state. Hence we consider that the effect of the pseudogap 
is crucial to understand the origin of the CO in hole-doped 
cuprates. In electron-doped cuprates, on the other hand, 
the CO tendency is observed in the paramagnetic state 0. 
This may be a reason why the present theory works for 
that case. 

While we have assumed the paramagnetic state, one 
may wonder whether the present theoretical framework 
actually predicts the strong asymmetry of the pseudogap 
between electron-doped cuprates {t' > 0) and hole-doped 
cuprates {t' < 0) simply by taking a different sign of t'. 
The origin of the pseudogap remains controversial even 
in hole-doped cuprates [351 ES] . If one assumes that the 

factor of 1/2 here comes from a large-formalism, where 
t is scaled by 1/A^. We may then invoke N = 2 when making a 
comparison with experiment. 

^Although a pseudogap was reported in the optical conductiv¬ 
ity spectra in the non-superconducting crystals m, the pseudogap 
corresponding to the observed one in hole-doped cuprates, namely 
a gap-like feature above the superconducting phase, is missing or at 
least very weak. 


pseudogap is driven by a strong charge-order tendency as 
hinted in some experiments [571|3i IMlIin], a theoretical 
study m using the present large-V scheme indeed sug¬ 
gests that pseudogap features should appear much weaker 
in electron-doped cuprates than hole-doped cuprates in 
the sense that a charge order tendency becomes much 
weaker in the former especially in a moderate doping re¬ 
gion relevant to the pseudogap (see fig. 6 in ED). 

5 Summary 

Motivated by the recent measurements by RXS and RIXS 
in electron-doped cuprates, we have studied charge exci¬ 
tation spectra associated with a d-wave bond order in the 
two-dimensional t-J model on a square lattice. We find 
that the static d-wave bond-order susceptibility Xd(q, 0) 
has two peaks at q, = (0.847r, 0.847r) and q2 = (0.497r, 0), 
which are generated by the 2kF scattering processes en¬ 
hanced by the d-wave character of bond-charge density. 
In spite of the proximity to the d-wave CO instability at 
q Ri qi, the peak at qi is very broad and becomes sharp 
only in the vicinity of its instability. On the other hand, 
the peak at q 2 becomes sharper with decreasing tempera¬ 
ture but does not diverge, indicating that the CO with mo¬ 
mentum q 2 is short ranged. These features are seen also 
in the equal-time correlation function S'(q). The spectral 
function of the d-wave bond order (Imxd(q, w)) forms a 
V-shape dispersion near q = (0, 0). This dispersion comes 
mainly from particle-hole excitations. The spectra bend 
back and reach close to zero energy at q = qi and q2 
where both the static d-wave susceptibility and the equal¬ 
time correlation function show a peak. The resulting spec¬ 
tra have charge gap-like features with a maximal gap at 
q Ri iqi and iq2. We argue that the CO observed in 
RXS is interpreted as a short-range order, which may not 
develop to become long range. It is interesting to explore 
gap-like features of the energy-resolved spectra and a pos¬ 
sible CO near q r; qi by RIXS and usual x-ray diffraction 
measurements, respectively, for electron-doped cuprates. 
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